Abstract. Molecular hydrogen is quite underabundant 
in damped Lyman-a systems at high redshift, when com- 
pared to the interstellar medium near the Sun. This has 
been interpreted as implying that the gas in damped 
Lyman-a systems is warm, like the nearby neutral inter- 
cloud medium, rather than cool, as in the clouds which 
give rise to most H I absorption in the Milky Way. Other 
lines of evidence suggest that the gas in damped Lyman-a 
systems - in whole or part - is actually cool; spectroscopy 

■ of neutral and ionized carbon, discussed here, shows that 
i: the damped Lyman-a systems observed at lower redshift 

z < 2.3 are largely cool, while those seen at z > 2.8 are 

■ warm (though not devoid of H2). To interpret the obser- 
5_( ' vations of carbon and hydrogen we constructed detailed 
Qh, numerical models of H2 formation under the conditions of 

' two-phase thermal equilibrium, like those which account 
, for conditions near the Sun, but with varying metallic- 
^f) ' ity, dust-gas ratio, etc. We find that the low metallicity of 
damped Lyman-a systems is enough to suppress H2 for- 
mation by many orders of magnitude even in cool diffuse 
clouds, as long as the ambient optical/uv radiation field 
is not too small. For very low metallicity and under the 
, most diffuse conditions, H2 formation will be dominated 
' by slow gas-phase processes not involving grains, and a 
; minimum molecular fraction in the range 10^^ — 10^^ is 
' expected. 
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1. Introduction. 

Perhaps it is something of a paradox, but large ground- 
based optical telescopes like Keck and VLT routinely do 
absorption-line spectroscopy of the gas in high-redshift ob- 
jects which surpasses what can be done for the interstellar 
medium (ISM) of our own Galaxy near the Sun. Viewing 
neutral gas at high redshift makes the dominant ion stages 
(C II, Fe II, etc.) accessible, to say nothing of Lyman se- 
ries in H I and D I and the Lyman and Werner bands of 
H2. The high spectral resolution and large collecting areas 
of modern ground-based instrumentation cannot easily be 



a low but non-negligible level 0.1 - 0.01 Solar, makes it 
possible (or, perhaps, merely hopeful?) to talk about the 
"interstellar medium" in these systems. The low metal- 
licity and general underabundance of dust (which is not 
directly obse rved in any one object but can be inferred 



statistically ( |Fall fc Pei, 1993| ; |Pei et al., 1991]) or from 



patterns of gaseous abundances ( Boisse et al., 1998 )) may 
render the gaseous medium in damped Lyman-a systems 
only more extreme versions of those in local dwarf systems 
like the LMC and SMC [| . But it has also been argued 
that the apparent recognizability of patterns in the ab- 



matchtjd in space-based iusUumeuLs, and ceiLaiuly nuL fui 
the 



sorption spectra has been over-interpreted (Izotov et al 



I same cusL 

Thus we may have exquisitely detailed and sensi- 
tive spectra of systems which cannot be imaged, whose 
nature 
terns c 



2001), and that the similarity of intermediate and low ion 



is Lheiefuie ML Lu be lufeiied hum Llieii paL- 



kinematics (Al III and C II or Fe II; C IV and C II behave 
differently) could mean that substantial ionization correc- 
tions to the metallicity are needed (however, see Vladilc 



Haehneh et al, 1998: Ledoux et al., 1998; McDon 



f gas kiuemaLlcs (Piucliaska Wolfe, 1997rl^§g 



aid fc Miralda-Escude, 1999) and ionization (Wolfe fc 



et al. (2001) for a contrary opinion). 



Prochaska, 2000a , b|) . So it is with most damped Lyman- 



a 

systems, defined as those absorption-line systems having 
N(II I) > 2 X 10^*^ cm"^ when seen against the emis- 
sion of background QSO's. Although well-formed galax- 
ian systems can and do harbor some of them at low red- 
shift, d amped Liyman a eystoms are for the most part bo — 
lieved |o bo protogalactic objects, porhape in disk eyfAoTue 



At the present time there seem to be several lines of 
evidence suggesting that, unlike the local ISM (where low- 
altitude neutral gas is perhaps 2/3 cool and 1/3 warm 
and the overall ratio including high-altitude material is 
1/2 and 1/2) the gas in damped Lyman-a systems is 
more predominantly warm. Examples include the high H I 
spin temperatures inferred from comparison of A21cm and 



Lyman-a absorption (Wolfe & Davis, 1979; Carilli ct al 



1996a, ^Chengalur & Kanekar, 2000; Kanekar & Chen- 



" ""J---", galur, 20011) (but see Lane et al. (2000)), the similarity of 
(Prochl aska fc Wolio) or perhaps m the ongomg merger of — 2 ' "1 ^ ^ F — — ^ y" . j • 1 



protogalactic clumps (Haehnelt et al., 1998). The numbers 
are such that damped Lyman-a systems seem to contain 
at least as many baryons as can be found in the local Uni- 
verse now, and it is of great interest to understand how 
these baryons become recognizable nearby objects within 
a relatively short redshift interval. 

The absorption-line gas in damped Lyman-a systems 
is easily discussed in the same terms, using the same phys- 
ical processes, that are employed locally. Evidence that 
the resi t - framc optical/uv radiation field is comparable to 



low and intermediate ion kinematics mentioned in the pre- 
ceding paragraph, and the very low column densities of II2 
discussed here. Along lines of sight with reddening EB_y> 
0.05 mag (N(II) > 3 x lO^*' cm^^) in local Copernicus 
spectra, it is always the case that a few percent or more of 
the neutral hydrogen is molecular. By contrast, the molec- 
ular hydrogen fraction in damped Lyman-a systems is 2-4 
orders of magnitude lower (i.e. 10"''' — 10~^), which can 
be interpreted as meaning that the gas temperature must 
be ab ove 3000 K ( Petitjean et al., 2000| ) Lanzetta et al 



that n ear the Sun (Lcvshakov et al., 2002 ;- Molaro ct al 



2002 ; Petitjean et al., 2000| ) and the presence of metals at 



(19891 ) used N(H2)/N(H) to constrain the dust/gas ratio 



Send offprint requests to: H. S. Liszt 
Correspondence to: hliszt@nrao.edu 



Boisse et al. (1998| ) show that there there is a bias against 



finding high metallicity damped Lyman-a systems at high N(H 
I), because the background sources are extinguished. 
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toward Q1337+113, similar to the approach taken in this 
work. 

Here we consider the formation of molecular hydro- 
gen in a gas which is in two-phase thermal equilibrium 
at low metallicity. Perhaps because H2 has been observ- 
able so rarely in the local ISM there is not a big literature 
on this subject, but the extant 1970's-era Copernicus ob- 



servations are well-explained in this way (Liszt & Lucas 
2000 |)T|ising modern sh ielding factors for radiative dissoci- 
ation (|Lcc ct al., 1996 ). The only surprise (if there indeed 



is one) is the low densities that are required to start H2 
formation locally, and the fact that even a "standard" H 
I cloud dSpitzer, 1978 ) should have a molecular fraction 
of 10-30% deep inside. Local diffuse clouds also have sur- 
prisingly high abundances of complex polyatomic species 



which follow immediately upon the presence of H2 (Liszt 

& LuCc,s, 1996; Lucas & Liszt, 1996), a phenomenon which gas at the Solar radius n, 



is not well understood but which can be used to account 
for the observed abundancess of simpler species such as 
CO dLiszt fc Lucas, 200^ ). 

In the course of this work, we consider the ionization 
and fine-structure excitation of carbon, and the inferences 
which may be drawn from observations of carbon in the 
damped Lyman-a systems. While the importance of car- 
bon to the physical state of the gas and the conditions for 
forming II2 cannot be stressed too highly, the extended 
discussion here in fact arose because of something of a 
coincidence. Searching the literature, it quickly became 
evident that there is a high degree of overlap between the 
two relatively scant datasets for carbon (chiefly, C I, C II, 
and C II*) and H2. 

The plan of the current discussion is as follows. In Sect. 
2 we lay out the basics of a calculation of two-ph ase equi- 
librium essentially following Wolfire et al. (1995a ) . We ex- 
amine the sensitivity of two-phase equilibrium to varia- 
tions in abundance, depletion, incident radiation and the 
like, in order to extract from the calculations those as- 
pects which can be related to existing absorption line data 
on damped Lyman-a systems. In Sect. 3 we compare the 
results of these calculations with observed hydrogen and 
carbon column densities: we show that the observations 
are consistent with an origin in largely cool gas for the 
damped Lyman-a systems at lower z (z < 2.3) and in 
warm gas at higher redshift. In Sect. 4 we describe calcu- 
lations of the H2 abundance in both warm and cool diffuse 
gas, under conditions of varying metallicty, etc., employ- 
ing (slow) gas-phase processes to explore the minimum 
expected amounts of molecular gas, and the more usual 
grain-catalysis (as in Liszt fc Lucas (2000| )) for cooler re- 
gions of higher molecular fraction. We also compute the 
variation of the molecular hydrogen fraction in small gas 
clots of constant density. From this, it follows that the low 
metallicities of damped Lyman-a systems are by them- 
selves sufficient to cause decreases in the molecular frac- 
tion by many orders of magnitude, even if cool neutral 
clouds are present. 



2. Two-phase equilibrium of the ISM 

2.1. Particulars and parameter sensitivities 

The ISM is so complicated that any model of it is bound 
to be heavily idealized and stylized, especially the assump- 
tion of a strict equilibrium. But some aspects, especially 
the presence of phases - discrete regimes of density and 
ionization - seem robust to variations in the underlying 
parameters. The existence of multiphase gas seems to be 
quite general, highly conserved across space and time. 

The basis of the present work is a calculation of two- 
phase equilibrium following work on the local ISM by 
Wolfire et al. (1995a). Notable constraints on the model lo- 
cally are the thermal pressure range P/k = 10^ — 10'' cm~'^ 
K observed {via C I and C I*) in neutral gas locally by 
Jenkins et al. (1983), and the electron density in warm 

0.02 cm-'^ ([Taylor fc Cordes 



1993). As noted by Wolfire et al. (1995a), strict two-phase 



equilibrium is an idealized goal toward which the ISM may 
tend, but true equilibrium is hard to attain, especially in 
warm gas where the cooling timescales are long. 

Basically, this "standard" model of the local diffuse 
ISM is driven by the heating due (in increasing order of 
importance in warm gas) to cosmic-rays, soft X-rays, and 
- the dominant mechanism in both war m and cool gas 



the p hotoelectric effect on small grains ( Bakes fc Tielens 



1994 ). The charge state of the small grains, their heating 



rates and the recombination of ions on them (which dom- 
inates over gas-phase processes at moderate densities in 
diffuse gas) are all very sensitive to the electron density, 
so local thermal and ionization equilbria are tightly cou- 
pled. Given that carbon is the main source of cooling in 
cooler gas (O I to a lesser extent) and an important source 
of free electrons, it follows that the equilibrium conditions 
and the condition of carbon in the gas are very tightly 
coupled. 



We recently used this code ( |Liszt, 200 1| ) to discuss the 
behaviour of the H I spin temperature in local warm gas 
(the H I there is not generally thermalized by collisions in 
multiphase equilbrium) and to calculate the abundances 
of H2 and CO in denser diffuse gas (Liszt fc Lucas, 2000). 



The latter reference describes the model in more detail 
than will be given here, especially the various coUisional 
processes which drive fine-structure cooling in the gas 
(they are, of course also discussed in the original refer- 
ence by Wolfire et al. (1995a)). 



Our basic calculation differs from that of Wolfire et al 



(1995a) in two minor ways which are subsumed by the ex- 



tent of our parameter variations. The earlier authors took 
the locally-determined soft X-ray flux from parmire et al 



(1992) and decomposed it into three components, each 



of which they represented by a (physically-motivated) 
plasma emissivity which was then separately attenuated 
(or not) by an assumed column of warm gas; the three 
contributions were then summed to provide the assumed 
incident X-ray spectrum at a typical location. Their final 
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Fig. 1. Ionization and thermal equilibrium calculations in atomic gas. The pressure P/k is shown as a function of the 
density of H-nuclei n(H), for four sets of parametric variations, a) Upper left; the metallicity, (dust/gas, C/H, 0/H 
etc.) varies in steps of 2 from 4 times to 1/32 times its reference value; b) Upper right, carbon and oxygen are removed 
(depleted) from the gas phase in steps of 2; c) Lower left, the "interstellar" (ambient) radiation field (ISRF) varies 
from 4 times to 1/32 times the reference value; d) Lower right, the impinging flux of soft X-rays is scaled. 



spectrum is characterized by a quantity they called Nw, 
the overall column of warm gas attenuating part of the 
incident spectrum ^ whose standard value they took to be 
Nw = 10^^ H atoms cm~^. We used a simpler representa- 
tion of the soft X-ray flux, whereby the observed spectrum 
is attenuated by a column of neutral gas Nw directly. Our 
standard value of Nw = 2 x 10^^ H atoms cm~^ produces 
a pressure-density curve which differs little from that of 



Wolfirc ct al. (1995a). 



We also used a somewhat different set of reference 
(Solar) atomic abundances whereby log(C/H) — -3.44, 
log(0/C) — 0.32, namely those which accompanied the 
distribution of the soft X-ray absorption cross- sections 
from ( Balucinska-Church fc McCammon, 199^ ), which 



both we and Wolfire et al. (1995a ) used (in their most 



^ The soft X-ray flux is so fragile that any cool gas interposed 
would absorb all of it, given that cool gas does not appear in 
small increments 



recently updated version). In both calculations the ref- 
erence model has (following [Wolfire ct al. (1995a )) little 
or no depletion of the gas phase oxygen and carbon onto 
grains, which might be acceptable for warm gas (where 
carbon does not bear the brunt of the cooling or con- 
tribute many of the electrons, see below) but is unlikely 
to be reasonable in local cool gas; if there are grains, they 
have to be made of something! The gas phase depletion is 
considered an adjustable parameter here, independent of 
the metallicity. 

The net effect of these two differences is slight. With 
the exception of one system, the observations show that 
carbon is present in the gas in the same proportion as 
other species which are typically undepleted locally, i.e. 
there is low metallicity, but little global carbon depletion. 

In Fig. 1 we show the basics of two-phase equilibrium 
and its sensitivity to some assumed parameter variations. 
In these curves of pressure vs. density (one assumes a 
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Fig. 2. As in Fig. 1, but the quantities shown in each panel are the energy loss rate nA due to radiation in the C II* 
A158/X cooling transition (lower, full curves) F, the total energy input to the gas (upper, shaded curves), both in units 
of erg . In cool gas (to the right), C II is by far the dominant coolant; in warm gas, most of the cooling occurs 



in the Ly-a hue. The mean C II cooling rate for the Milky Way is 2.65 x 10 ° erg s H ^ (Bennett ct al., 1994) 



density, calculates the ionization and thermal equilibrium, 
sums the particle densities and multiplies by the derived 
temperature), the regions where dP/dn < are held to be 
unstable and unlikely to occur in nature. Those pressures 
for which there are two densities n' having dP/dn|„' > 
are those for which "two-phase" equilibrium is possible. 
Two-phase heating and cooling calculations by themselves 
furnish only the possib ility of multi-phase equili brium; 
Hennebelle (2000| ) and [Kritsuk fc Norman (200^) elab- 
orate on some processes by which phase transitions and 
multi- phase equilibrium are actually brought about. 

The curve labelled "1" in all panels is the same (ref- 
erence) model; perhaps most clearly at upper right it is 
apparent that the reference model provides for two-phase 
equilibrium over precisely the pressure range which is ob- 



served in local gas (Jenkins et al., 1983). Typically (but 



see below and the lower- left panel), only warm (7000 K 
- 10000 K) neutral gas appears if the density and pres- 
sure are below the two-phase regime and only cold (below 



1000 K) neutral gas appears if the pressure and density 
are larger. If conditions are such that the warm gas may 
persist up to higher density and pressure, or when the 
cold medium cannot exist except at higher density and 
pressure, the medium is more likely to be warm. 

At the top left, the metallicity is varied by factors of 4 
to 1 /32 relative to the reference value ; these may be com- 



pared with the calculations of Wolfire et al. (1995b) . Many 
properties of the model actually vary in concert including; 
the elemental and gas-phase C and O abundances (there is 
no depletion); the columns of all other species heavier than 
Helium providing soft X-ray absorption; and the number 
of grains both large and small (the latter scales the heating 
rate of photoelectric ejection). Vladilo (1998) shows that 
the dust/metal ratio D/Z « 0.6 at higher redshift rela- 
tive to conditions in the Milky Way (MW) so that D/G = 
D/Z X Z/G changes mostly because of metallicity (Z/G), 
not D/Z. In fact, such an overall scaling has a surprisingly 
modest effect compared to some other changes discussed 
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Fig. 3. As in Fig. 1 and Fig 2, but the quantity shown in each panel here is the ratio of once-ionized to neutral carbon 
n(C II)/n(C I). 



next. The warm phase may persist up to slightly higher 
pressures and densities in gas of lower overall metallicity 
but the existence of two-phase equilibrium is not threat- 
ened as long as the balance between heating (small grains) 
and cooling (carbon) is maintained in denser gas (as it is 
in this case). 

At upper right, the gas-phase abundances of C and O 
are decreased while all other parameters are held fixed and 
this actually has a much more profound effect than scal- 
ing the grain (heating) and gas (cooling) abundances in 
tandem. Carbon provides some heating of the gas {via its 
ionization, see Spitzer fc Scott (1969| )) but is more impor- 
tant as a source of electrons, influencing the grain heating 
and recombination rates, and as a coolant. Clearly, when 
the ambient fluxes are maintained while the main coolants 
(carbon and oxygen) and source of electrons are depleted, 
it is substantially more difficult to provide enough pres- 
sure to maintain gas in the cool phase. This hearkens back 
to our discussion, slightly earlier, of the lack of inclusion 
of depletion in the reference model for MW conditions. In 
order to maintain two-phase stability over the reference 



range of pressures in the presence of gas-phase depletion, 
the incident optical/uv ionizing radiation field must be de- 
creased somewhat. In fact this might not be unreasonable 
locally for cool gas, which does not occur in arbitrarily 
small increments and therefore exists under conditions of 
non-trivial extinction (especially at uv wavelengths). At 
lower metallicity, this sort of compensation is less obvious 
and (consistent with some aspects of the observations) the 
pressure may in fact be somewhat higher in cool gas. 

At lower left, we scale the ambient optical/uv radiation 
field which ionizes carbon (directly) and provides heating 
and electrons via the photoelectric effect on small grains. 
Regions whose ambient radiation field (ISRF) falls even a 
factor four below the average will be predominantly cool. 
This is an inversion of the normal order of things in diffuse 
gas and darker clouds in the nearby ISM occur at higher, 
not lower, thermal pressure. 

Last, at lower right, we scale the soft X-ray flux. Cool 
regions of higher density will probably never see the full 
extent of the ambient soft x-ray flux because even a small 
region of appreciable density accumulates a high degree 
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Fig. 4. Fine-structure excitation calculations for varying 
Tcmb= z' X 2.73 K. Top: the ratio of first-excited to ground- 
state populations in neutral carbon, CI. Middle, the ratio 
of excited to ground state populations in once-ionized car- 
bon, C II. Bottom, as in the middle panel but for the case 
where carbon is depleted by a factor 4 from the gas-phase, 
relative to the standard model. 



of X-ray absorption. The heating due to soft X-rays tends 
to increase on a per-event basis as the X-ray spectrum 
hardens somewhat after absorption, and tends to decrease 
as the electrons cause more secondary ionization and less 
direct heating when the ionization fraction is smaller (in 
neutral gas). 



2.2. Thermodynamics 

Fig. 2 shows, for the same set of parameter variations, the 
total gas heating rate F and the loss of energy (nA) due 
to cooling in the A158/i C II line (discussed below at some 
length), in units of ergs s~^ H^^. For the reference model, 
some 90% of the heating is due to the photoelectron ef- 



fect operating on small grains, see Fig. 3 of Wolfire ct al 



(1995a). In warm gas, the burden of cooling is taken up 
by excitation of the Ly-a lines under the highly-idealized 
assumption that all photons escape; in fact, these pho- 
tons may travel substantial distances but by and large do 
not always escape a two-phase galactic layer before be- 
ing absorbed by dust, at Solar metallicity: the situation 
may actually b e ver y different i n three-phase models, see 
Neufeld (1991 ) and Liszt (2001 ). In cool gas, the cooling 
is due almost entirely to carbon and oxygen fine-structure 
excitation, in the reference model. The curves of nA in 
Fig. 2 directly show the brightness (normalized; per H) 
of the C II* A158/X cooling transition ; in cool gas the O 
I* transition at A63/i makes up most of the cooling not 
provided by C II. 

There is really a quite profound change in the lo- 
cal thermodynamics when the metallicity or ionizing flux 
varies (Fig. 2, upper and lower left). In particular, while 
the energy input per H into cool and warm gas is nearly 
the same for the "Solar" metallicity and the standard 
ISRF or higher, it is much smaller in cooler gas when 
the metallicity is low. In the Milky Way, the carbon cool- 
ing rate can be used to infer the heating in both cool and 
warm gas. In systems of low metallicity such is not the 
case. 

In general, the brightness of the C II 158/i line scales 
with metallicity at all densities (Fig. 2, upper left) and 
varies nearly linearly in cool phase-stable gas since n(C 
II*)/n(C II) changes little with metallicity (see Fig 5.). 
It also scales with the strength of the ISRF in cool gas 
(Fig. 2, lower left), due to changes in n(C II*)/n(C II) 
(Fig. 5). The C II A158/i brightness varies with depletion 
at lower densities (Fig. 2 upper right) and slightly, at low 
density, with variations in the soft X-ray flux. From the 
results at upper right in Fig. 2 we sec that C II may not be 
the dominant host of cooling in cooler gas, at quite high 
densities, if the depletion (not the metallicity) is extreme; 
oxygen becomes the preferred coolant due to the higher 
energy separation in its ground state flne-structure levels. 
There is no reason to believe that such conditions occur 
widely in diffuse gas, either locally or in damped Lyman-a 
systems, although one source in Table 2 (PIIL957) seems 
quite deficient in carbon, given its quoted metallicity. 

2.3. Carbon ionization 

Fig. 3 shows the ratio of the two lowest ionization states 
of carbon n(C II)/n(C I). It differs by a factor of 40 - 
100 between warm and cool neutral gas for the standard 
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model, and is relatively insensitive to parameter changes, 
redshift, etc. This makes it a sensitive indicator of the 
thermodynamic state of the g discussed in Sect. 3, 

but also causes confusion when mixtures of the two phases 
are observed along the same line of sight (see Sect 2.6). 
Some aspects of the behaviour shown in Fig. 3 are counter- 
intuitive; for instance, a stronger ISRF maintains a high 
C II/C I ratio into denser gas, but the C II/C I ratio is 
actually smaller in very tenuous gas when the radiation 
field is higher. This arises in part because the equilibrium 
pressure increases with the strength of the radiation field 
and partly due to the level of ionization in the gas (which 
is not solely determined by carbon). 

If the C II lines are saturated or N(C II) is other- 
wise unknown, but N(C I) and N(H) are available, the 
thermodynamic state of the gas may also be inferred by 
computing N(C)/N(H) Z N(C II)/N(H) « N(C II)/N(C 
I) X N(C I)/N(H). Presumably, if the gas is actually cool, 
the inferred carbon abundance will be much too large (in- 
consistent with the metallicity) if the high n(C II) /n(C I) 
ratios typical of warm gas are assumed. 

2.4- C I fine- structure excitation 

The triplet fine-structure of the C I ground state is also 
accessible to spectroscopy. It is actually rather insensitive 
to the parameter variations noted above (with the slight 
exception that C I*/C I increases modestly with increasing 
X-ray flux at low density) but responds very strongly to 
a change in the cosmic background level because the first 
excited state is only 23.2 K above ground. Fig. 4 (top) 
shows the ratio of populations in the lowest two levels as a 
function of Tcmb- Clearly, this ratio is a sensitive indicator 
of the ionization state at low redshift (it differs by a factor 
of about thirty between warm and cool gas) but is rather 
insensitive for redshifts above 1.5 or so. The C I*/C I 
ratio measures Tcmb at high redshift in the context of 
our models (see Roth & Meyer (1992)). Although perhaps 
within observational errors, it follows from the results in 
Fig. 4 that quoted column densities of neutral carbon for 
high-z systems should include a correction factor for the 
excited state population. 



2.5. Cooling rates and C II fine- structure excitation 

The rate of energy loss due to a radiative transition be- 
tween two levels k and j having statistical weights gk and 
gj, level populations n^ and nj (units of cm^"^), energy 
separation E^j and spontaneous emission rate Akj (s"^), 
immersed in an isotropic radiation field characterized by 
Tcmb is dGoldreich fc Kwan, 197^ 



nAkj = ^kj/3kj-Ekj(nk - ("jgk/gj - nk)Fkj (Tcmb)) 

where Fkj(x) = l/(exp(Ekj/kx)-l). /3kj is the probabil- 
ity of photon escape, typically varying as /3kj = (l-exp(- 



UTkj))/(uTkj) where u is a geometrical factor of order unity 
and Tkj is the optical depth of the transition. The escape 
probabability is taken as unity in these two-phase mod- 
els, which is acceptable because the various fine-structure 
transitions do not achieve high optical depth in typical 
diffuse clouds. In calculation of the local statistical equi- 
librium, these expresssions allow for a proper evaluation 
of the fine-structure cooling rates for varying Tcmb- 

The 133.5 and 133.6 nm C II and C II* absorption 
lines have received more attention because of the obvious 
possibility that C II is the dominant ionization stage. But 
discussions of the excitation of the C II* fine-structure 
levels often fail to recognize the role of carbon in deter- 
mining the conditions under which it is observed; in many 
cases (see Fig. 2) C II will be the dominant coolant when 
it is the dominant ionization stage. This means that the 
brightness of the C II A158Ai fine and the C II*/C II ra- 
tio are determined by thermodynamic equilibrium - the 
energy input to the gas which occurs at a given density 
and temperature. So, if the energy input into the gas does 
not change, neither will the brightness of the A158/Lt line 
or the strength of the 133.5 nm C II* absorption line. In 
Fig. 2 at upper right, for a density of n(H) = 20 cm~'^, 
it requires a factor 8 drop in the the amount of carbon in 
the gas to weaken the A158/i line by about a factor 2, all 
other things being equal. 

Fig. 5 shows the parameter sensitivities of the C II*/C 
II ratio. Weak variation in the n(C II*)/n(C II) ratio in 
cool gas is shown in many ways in the panels of Fig. 5. For 
the reference model, the ratio varies by only 50% between 
n(H) = 2 cm-3 (Tk= 1000 K) and n(H) 40 cm'^ (Tk= 75 
K). The variation with carbon depletion at upper right il- 
lustrates the effects of thermodynamic considerations. The 
ratio n(C II*)/n(C II) increases as carbon is removed from 
the gas or n(C II) decreases, because the brightness of the 
A158/i line (and n(C II*)) must stay constant in order to 
carry away a given amount of energy which is being in- 
put to the gas (see Fig. 2). The behaviour at lower left, 
where the C II* /C II ratio is a scaled representation of 
the cooling rate (also see Fig. 2) shows how the brightness 
of the A158/Lt line, or the amount of excited-state C II*, 
must change proportionally as the energy input to the gas 
changes. 

In warm gas where its excitation is weak and C II is 
not the dominant coolant, the population of the C II* 
level is susceptible to change with Tcmb- This is shown in 
the lower two panels of Fig. 4. The ability of the C II* 
population to discriminate between phases may be lost at 
higher redshift. 

In this context it is important to note that the AE/k 
= 91 K energy separation of the C II fine-structure levels 
(to a lesser extent the 232 K separation in O I), is ab- 
solutely crucial to the similarity of two-phase equilibrium 
conditions out to moderate redshift, say z < 10. If the 
gas locally were cooled by softer photons of energy AE 
such that we could observe out to some redshift where 
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AE/kTcmb^ 4, the distant gas would have to have a dif- 
ferent coolant. It is apparent from the behaviour of the C 
I*/C I ratio in Fig 4 that C I*, for example, could not be 
an important gas coolant at z > 3. The same is true of 
the lower lines of CO (where AE/k = 5.5 K, 11.0 K, 16.5 
K, etc.). 

2.6. Observing admixtures of the two phases 

When a property like the N(C II*)/N(C II) ratio - caU it 
R - takes on values Rw and Rc in warm and cool gas, 
respectively, and a transparent superposition of phases is 
observed to have a global value Robs, Robs is related to 
the proportion of gas in the two phases as 



Fc 



Rr 



Rw 



Rw — Rc Robs 
where Fc is defined as the fraction of the gas in the 
cool phase. 

In stable two-phase gas (Fig. 3 at upper left) Rw= 
20000 and (roughly) Rc = 700/n(H) in gas of Solar metal- 
licity (Rc varies rather less at low metallicity) . In this 



case, measuring a global average of, say, Robs — 1000 im- 
plies only that there is some fraction Fc < 1 of gas having 
Rc< 1000 along the line of sight. 

This ambiguity cannot generally be resolved with re- 
course to comparison of profiles, which are largely undif- 
ferentiated in the C II line, having been broadened due to 
some combination of; underlying kinematics in the host 
system; doppler-broadening of individual gas complexes; 
and, saturation. In the Milky Way, C II lines are typically 
fit with much larger b- values (Doppler broadening param- 
eters) even when a single cloud or ga s complex dom inates 
the line of sight to a nearby star (c/. Morton (1975 )). Be- 
cause of this one cannot, for instance, try to compare the 
C II* with only that fraction of the C II which overlaps it 
in velocity. 

3. Warm or cool gas in damped Lyman-a 
systems? observations of C I and C II 

In Tables 1 and 2 we have gathered information on the 
neutral and ionized carbon lines available in the litera- 
ture. For each source we show the observed column den- 



10 H. Liszt: H2 formation in damped Lyman-a systems 

Table 1. Hydrogen and carbon at lower redshift 



Source 


3C286 


0454+039 


1756+237 


1331+170 


0013-004 


1232+0815 


PHL957 


z 


0.69 


0.86 


1.67 


1.78 


1.97 


2.24 


2.31 


[Zn/H] 


-1.22 


-1.1 


(-0.8)" 


-1.27 


-0.80 


-1.20 


-1.38 


N(H I) 


1.8E21 


4.2E20 


2.0E20 


1.5E21 


5.0E20 


7.9E20 


2.8E21 


N(H2) 










6.9E19 


6.3E16 


<5.0E15 


N(C I) 


>4E13 


4.4E13 


1.4E13 


1.7E13 


6.3E13 


4.0E13 




N(C II) 








4.5E16 


2.7E16 




7.8E15 


N(C II*) 






7.9E13 


1.4E14 


1.9E14 




7.6E13 


N(C 11)^" 


3.8E16 


1.2E16 


(1.1E16) 


2.9E16 


2.8E16 


1.8E16 


4.1E16 


N(C II)/N(C 1)^= 


(<952) 


(293) 


(803) 


2647 


446 


445 


>1734 


N(C II*)/N(C II) 






(0.0072) 


0.0031 


0.0070 




0.0097 



assumed for illustrative purposes; other (entries) follow 
'' N(C II)m is the amount of carbon if [C/Zn] is Solar 

values in parentheses use N(C II)m 
References: 
3C286: 



Boisse et al. 



Boisse et al. (199^); iRoth fc Meyer (19921) 



(19981): 



0454+039: 
1756+237: 

1331+170: IChaffee ct al. (19881): |Ge et al. (1997D 



Roth fc Bauer (199d ) 



0013-004: |Ge et al. (20011) 



Steidel et al. (1995) 



1232+08 15: [Black et al. (19871) ; [Srianand et al. (2000| ) quote N(H2) 
PHL957: [Black et al. (1987|) 



1.5E17 cm" 



sities, and the metallicity derived from zinc, which is be- 



lieved to be (at most) only lightly depleted (Pettini et al 
1994; ^ladilo et al., 2000). From the metallicity we derive 



a quantity N(C II)m which is meant to represent the max- 
imum amount of carbon available in the gas phase, i.e. it 
is the column density of carbon assuming that carbon and 
zinc have the same metallicity and depletion. Comparison 
with the measured C II column densities shows that car- 
bon is not lacking, with the notable exception of PHL957. 
N(C II)m is used as a surrogate where measured values of 
N(C II) are unavailable at low redshift in Table 1. 

The ratio N(C II) /N(C I) is much larger for the sources 
at high z in Table 2, above 13,000-30,000, consistent only 
with an origin in overwhelmingly warm gas: the fraction of 
cool gas must be of order a few percent at most. At lower 
redshift, two of the measured ratios, and that inferred for 
the source 1756+237, are below 1000, consistent with cool 
neutral gas in the stable region of the two-phase equilib- 
rium, at thermal pressures P/k = 3,000 - 4000 cm~'^ K, 
and densities n(H) = 10-70 cm^^. The data for PHL957 
are also consistent with stable cool gas at n(H) ^ 10 cm^^, 
p/k 2500 cm-3 K, while 1331+170 is in the marginally 
unstable region at n(H) = 10 cm~^, p/k = 2800 cm~^ 
K. For the systems at lower redshift, there must be a very 
large proportion of cool gas, though we cannot say exactly 
how much. 

The N(C II*)/N(C II) ratios are noticeably smaller at 
high redshift in Table 2, and would presumably be smaller 
still by a factor 2 or more, were it not for the cosmic 
background radiation, which has a noticeable effect (only) 
above z = 2.5 in warm gas (see Fig. 4). The N(C II*)/N(C 



II) ratios observed at lower redshift, in presumably cool 
gas, are 50-100% higher than in the models. This could 
be explained by an enhanced ISRF (Fig. 5, bottom left 
panel) or by some depletion of carbon for PIIL957, where 
the metallicity of carbon is (uniquely) very low compared 
to that of Zn, and N(C II*)/N(C II) is rather large. 

To summarize, comparison of the tabulated observa- 
tional results with the calculated properties of carbon in 
two-phase media shows quite unambiguously that the ob- 
served gas has a substantial contribution from cool gas at 
lower redshift but is very largely warm at z > 2.8. Whether 
this shift from cool to warm g as is systematic or coinci- 
dental remains to be explored. Norman fc Spaans (1997 ) 
predicted that multi-phase equilibrium, cool gas and sub- 
stantial quantities of H2 would occur in protogalactic disks 
only when the metallicity had increased to 0.03 - 0.1 Solar, 
somewhere in the interval 1 < z < 2. 

4. The formation of molecular hydrogen 

Levshakov et al. (2002) provide a summary of the molec- 
ular fractions observed toward 14 sources (their Table 2) 
including most of those discussed here. The lowest molecu- 
lar fraction is seen toward 0000-262 at z = 3.39, the high- 
est toward 0013-004 at z = 1.97 (see Tables 1-2 here). 
The latter is the only direction for which the molecular 
fraction exceeds 1/3000, and most are below 10^^. In this 
Section we address the existence of damped Lyman-a sys- 
tems with large hydrogen column densities, substantial 
amounts of cool gas, but very small fractions of molecular 
hydrogen. The models of Il2-formation presented here are 
very similar to those of Liszt fc Lucas (2000 ), who consid- 
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Table 2. Hydrogen and Carbon at highest redshift 



Source 


1337+113 


0528-250 


0347-382 


0000-262 


z 


2.80 


2.81 


3.03 


3.39 


[Zn/H] 


-1.00 


-0.91 


-1.23 


-2.07 


N(H I) 


8.0E20 


2.2E21 


2.52E20 


2.6E21 


N(H2) 


<5.0E16 


6.0E16 


8.2E14 


1.1E14 


N(C I) 


<1.6E13 


<5.9E12 


<4.0E11 




N(C II) 


2.0E17 


1.7E17 


5.1E15 




N(C II*) 




3.6E14 


1.9E13 




N(C II)n. 


2.8E17 


9.6E16 


3.1E15 


7.9E15 


C II/C I 


>12658 


>28862 


> 12700 




C II*/C II 




0.00211 


0.00389 




References: 
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III nil 1 1 1 1 nil 
ISRF 


1 1 nil 1 1 1 1 III 
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CM 
















^ 10-2 
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10-9 










1 


III mil 1 1 1 1 mil 


1 1 mil 1 1 1 1 iiir 



Lanzetta et al. (1989 ) 

19971); iSrianand 



1337-H13 
0528-250: pc ct 

0347-382: Lcvshakov et al. (2002|) 
0000-262: Prochaska & Wolfe (19991) 



Petitjean (1998) 



Levshakov et al. (2000) 



of early-universe conditions must be modified to include 
relevant values for the photodissociation of H~ and H2 
and the cosmic-ray ionization of H2. For the latter we as- 
sumed C = 2x~^'^s~-'^ per H and 1.08 C per H2), and for the 



txxrr 



trr 



T 



-ftr 



photodissociation rate of H2 in free space we followed Lee 



n( l I ) [cm--"*] 

Fig. 6. Molecular fractions 2n(H2)/n(H) in free space re- 
sulting solely from slow gas-phase H2-formation mecha- 
nisms, for the two-phase models and paramter variations 
shown at left in Fig. 1. See Sect. 4 of the text. 



ered the formation of CO in local diffuse gas, given certain 
other conditions like the presence of £[00+, but there are 
a few differences which we now remark. 



J^.l. }i2 -formation in the gas phase 

The present calculations include X-ray heating and ioniza- 
tion, because they are integral to the question of two-phase 
equilibrium. X-rays have little effect on the large neutral 
gas columns which harbor appreciable molecular column 
densities nearby in the Milky Way but are also included 
now because we are interested in understanding the mini- 
mum molecular fractions which can be expected, and these 
are presumably set by the slow gas-phase processes (not 
involving grains) which formed the first stars, and which 
occur in low-density regions of (unshielded) free space. 

Gas phase H2-formation occurs via the exothermic re- 
action pathways H -|- e~ ^ H~ + hi', H~ + H ^ H2 -f e~ 
and H+ + H ^ H+ -h hv, H+ + H ~> H2 + H+ . Many of 
th e basic reactions are c ited by puy et al. (1993 ) and most 
by [Haiman et al. (1996| ) (aU can straightforwardly be lo- 



et al. (1996D . In order to formulate a treatment of the vari- 
ation of the photodissociation of H~ with extinction, we 



integrated the cross-sections of Wishart (197£ ) over the lo- 
cal ISRF, finding an unshielded rate of 1.5 x 10^^ s^^, just 
over half of which arises at wavelengths beyond 800 nm; 
the photo-dissociation is, therefore, not strongly attenu- 
ated under the conditions discussed here and we elected 
to ignore extinction in this regard (the rate quoted in the 
UMIST database is 2.4 x 10"'' s~^ but we used the smaller 
value). 

Current values of the reaction rates for all impor- 
tant processes are given in the UMIST reaction database, 
whose values we employed unless otherwise noted. Species 
followed during modelling of the chemistry included H I, 
He"*", H"*" and e~ - all given by the calculations of two- 
phase equilibrium -as well as H~ , H2 and . 

Fig. 6 shows the free-space abundance of H2 arising 
solely from gas-phase processes in the two-phase models 
under conditions of varying metallicity and ISRF (the two 
left-hand panels of Fig. 1). Unlike grain formation scenar- 
ios, which take advantage of high N(H) to boost the molec- 
ular fraction at high density, free-space gas-phase forma- 
tion of H2 does not seem to distinguish between warm 
and cool conditions. The molecular fractions calculated in 
Fig. 6 correspond well with the smallest values in the lo- 
cal ISM, or, for that matter, in damped Lyman-a systems. 
Unfortunately, this minimum is not diagnostic of the host 
gas conditions. 



cated in the UMIST reaction database) but the discussion 
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Nh=4x1 



So 



cm 



Gr=0^ 




0.2 0.4 0.6 0.8 1.0 
fractional distance into cioud 

Fig. 7. Molecular fraction as a function of fractional dis- 
tance into cool clouds of density n(H) = 32 cm~^ and 
various central column density N(H) as indicated. The up- 
permost curve in each panel is for local conditions. The la- 
belled curves are for varying metallicity and a dust /metal 
ratio D/Z = 0.6 (60% of the local value). The bottom 
curve represents formation solely in the gas phase, in the 
limit of zero metallicity parameter; see Sect. 5 of the text. 



4-. 2. The H2 abundance in cool clouds 

H2 formation in cool neutral gas clouds is illustrated in 
Fig 7, which indeed shows why even damped Lyman-a 
systems with appreciable cool gas still may lack molecu- 
lar hydrogen. To create this diagram, we considered (fol- 



lowing Liszt & Lucas (2000)) a spherical clot of gas of 
constant density, immersed in isotropic radiation fields 
(X-ray, cosmic-ray, optical/uv, etc.). This was compu- 
tationally divided up into 128 radial zones, in each of 
which we derived the temperature/ionization structure 
and H2 abundance. The latter requires iteration, because 
the maintenance of H2 is a sharply non-linear process de- 
pendent on the H2 column and extinction between any 



point and free-space (Lee et al., 1996). We adopted a fairly 



straightforward relaxation method which converged with 
gratifying rapidity. 

Calculation of the abundance of molecular hydrogen 
is typically made feasible by employing a set of shielding 
factors which account in an average way for the many very 
complicated effects of line-overlap and radiation transport 
in the dissociation process. We used the shielding factors 
of ( |Lee et al., 199^ ) which were calculated for local gas. 
The justification for this is that the dominant effect requir- 
ing consideration here is the order of magnitude change 
in the number of grains at very low metallicity, not the 
factor of two differences in individual grain properties be- 
tween local grains and those seen, for in stance in t he Mag- 
ellanic clouds. The parametrization of Pei (199^ ) for the 
Milky Way, LMC and SMC shows that, for a given amount 
of B-band extinction, the grain distribution provides suc- 
cessively somewhat more extinction at (say) Ly-a as the 
metallicity declines; the inference is that graphite grains 
disappear and silicates do not. But this effect is dominated 
by the overall diminution of the extinction with lowered 
metallicity. 

Figure 7 shows the radial variation of the fraction of 
H-nuclei in molecular form over spherical gas clots of con- 
stant density n(H) = 32 cm~'^ for different column densi- 
ties N(H) through the center of the clot. The mean line of 
sight averaged over the circular face of such a body inter- 
sects it at an impact parameter of 2/3 of the radius (at a 
value 0.33 along the horizontal axis in Fig. 7), where the 
column density is 3/4 of that through the center. 

In each panel of the Figure there are 8 vertically- 
separated curves. At top, shaded, is the result which 
would apply in the Milky Way, where we have taken the 
dust/metal ratio as observed locally (the reference model 
of the two-phase calculations) and depleted carbon and 
oxygen in the gas phase by a factor of 2.4. The bottom 
curve, also shaded, is the result when the metallicity goes 
to zero and only the gas-phase formation of H2 is included; 
a modest amount of self-shielding occurs and the molec- 
ular fraction is slightly higher than in free space (Fig. 6). 
The intermediate curves assume a dust/gas ratio 0.6 rel- 
ative to the re ference model for the Milky Way following 
Vladilo (1998| ) (a small effect at higher N(H) but of real 



importance to the thinnest model) and no depletion of car- 
bon and oxygen. These curves are labelled by their varying 
metallicity as in the previous diagrams. 

The cloud with N(H) ~ 4 x 10^° cm~ ^ would be a com- 
pact (4 pc), cool (130 K) ^pitzcr (197^) "standard" cloud 
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in the Solar neighborhood. It would also be very substan- 
tially molecular if found in the Solar vicinity, but a factor 
4 decline in the metallicity suffices to reduce the molecular 
fraction by some four orders of magnitude. Even the model 
having a four times higher column density (compare with 
the entries in Tables 1-2) cannot sustain an appreciable 
molecular fraction when the metallicity is reduced by a 
factor 10, which is hardly extreme for one of the damped 
Lyman-a systems. 

4-3. The role of geometry 

The role of geometry can also be inferred from Fig. 7. At 
any given metallicity, a cloud with lower N(H) produces 
much less than one-fourth as much H2 as that illustrated 
in the next-lower panel. This is another reason why the 
molecular fraction may vary widely between two lines of 
sight with similar N(H), N(C II)/N(C I), and/or metalHc- 
ity (for example). Molecular hydrogen will readily popu- 
late a region when the circumstances are propitious, but 
can easily be prevented from forming by the vagaries of 
local source structure. 

44. PKS0528-250 

This source (Table 2) has an overall molecular fraction 
5.5 X 10^^ despite the lack of evidence (in carbon) for 
any appreciable amounts of cool gas; earlier we asserted 
that (very roughly) no more than a few percent of the 
gas could be cool. [Carilli et al. (1996b| ) did not detect 
21cm H I absorption, placing a 2 — cr upper limit N(H I) 
< 2 X 10^^ cm~^Tsp. So, the molecule-bearing gas must 
be cool, occupying roughly 1% of the total gas column for 
Tsjp ^ 100 K. In the context of our models the gas must 
also be fairly dense, n(H) ^ 100 cm~'^, occurring over only 
a very small fraction of the path length (1 kpc or more ) 
occupied by the g whole. 

5. Summary 

The thermodynamics of the ISM near the Sun are strongly 
influenced by the disposition of carbon, in grains and in 
the gas. Spectroscopy of neutral and ionized carbon affords 
the opportunity to probe the processes which are most 
basic to the structure of the gaseous medium, nearby and 
in damped Lyman-a systems. For whatever reason, the 
very distant gas is easily understood in the same terms 
as that seen nearby. Here, we have shown that one seem- 
ingly disparate aspect of high-z systems, their small frac- 
tions of molecular gas, can also be easily understood. Even 
when cool gas is present, which must be the case for the 
systems discussed here at z < 2.3, abundances of H2 are 
suppressed by many orders of magnitude at lower metal- 
licity as a result of the sharply non-linear nature of the 
processes required to maintain substantial columns of II2 . 
No wholesale reorganization of the gaseous medium need 



be hypothesized to account for low H2 abundances. Con- 
versely, we showed that the slow gas-phase processes which 
formed H2 in the early Universe provide for a minimum 
molecular fraction in the range 10^^ — 10~^. 
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Abstract. Molecular hydrogen is quite underabundant in damped Lyman-a systems at high redshift, when 
compared to the interstellar medium near the Sun. This has been interpreted as implying that the gas in damped 
Lyman-a systems is warm, like the nearby neutral intercloud medium, rather than cool, as in the clouds which give 
rise to most H 1 absorption in the Milky Way. Other lines of evidence suggest that the gas in damped Lyman-a 
systems - in whole or part - is actually cool; spectroscopy of neutral and ionized carbon, discussed here, shows 
that the damped Lyman-a systems observed at lower redshift z < 2.3 are largely cool, while those seen at z > 
2.8 are warm (though not devoid of H2). To interpret the observations of carbon and hydrogen we constructed 
detailed numerical models of H2 formation under the conditions of two-phase thermal equilibrium, like those 
which account for conditions near the Sun, but with varying metallicity, dust-gas ratio, etc. We find that the low 
metallicity of damped Lyman-a systems is enough to suppress H2 formation by many orders of magnitude even 
in cool diffuse clouds, as long as the ambient optical/uv radiation field is not too small. For very low metallicity 
and under the most diffuse conditions, H2 formation will be dominated by slow gas-phase processes not involving 
grains, and a minimum molecular fraction in the range 10"** - 10"^ is expected. 



Key words, quasars: absorption lines; ISM: molecules 

1. Introduction. 

Perhaps it is something of a paradox, but large ground- 
based optical telescopes like Keck and VLT routinely do 
absorption-line spectroscopy of the gas in high-redshift ob- 
jects which surpasses what can be done for the interstellar 
medium (ISM) of our own Galaxy near the Sun. Viewing 
neutral gas at high redshift makes the dominant ion stages 
(C II, Fe II, etc.) accessible, to say nothing of Lyman se- 
ries in H I and D I and the Lyman and Werner bands of 
II2. The high spectral resolution and large collecting areas 
of modern ground-based instrumentation cannot easily be 
matched in space-based instruments, and certainly not for 
the same cost. 

Thus we may have exquisitely detailed and sensitive 
spectra of systems which cannot be imaged, whose na- 
ture is therefore left to be inferred from their patterns of 
gas kinematics (Prochaska & Wolfe, 1997, 1998; Haehnelt 
et al., 1998; Ledoux et al., 1998; McDonald & Miralda- 
Escude, 1999) and ionization (Wolfe & Prochaska, 
2000a, b). So it is with most damped Lyman-a systems, 
defined as those absorption-line systems having N(II I) 
> 2 X 10^*^ cm~^ when seen against the emission of back- 
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ground QSO's. Although well-formed galaxian systems 
can and do harbor some of them at low redshift, damped 
Lyman-a systems are for the most part believed to be pro- 
togalactic objects, perhaps in disk systems (Prochaska & 
Wolfe) or perhaps in the ongoing merger of protogalac- 
tic clumps (Haehnelt et al., 1998). The numbers are such 
that damped Lyman-a systems seem to contain at least 
as many baryons as can be found in the local Universe 
now, and it is of great interest to understand how these 
baryons become recognizable nearby objects within a rel- 
atively short redshift interval. 

The absorption-line gas in damped Lyman-a systems 
is easily discussed in the same terms, using the same phys- 
ical processes, that are employed locally. Evidence that 
the rest-frame optical/uv adiation field is comparable to 
that near the Sun (Levshakov et al., 2002; Molaro et al., 
2002; Petitjean et al., 2000) and the presence of metals at 
a low but non-negligible level 0.1 - 0.01 Solar, makes it 
possible (or, perhaps, merely hopeful?) to talk about the 
"interstellar medium" in these systems. The low metal- 
licity and general underabundance of dust (which is not 
directly observed in any one object but can be inferred 
statistically (Fall & Pei, 1993; Pei et al., 1991) or from 
patterns of gaseous abundances (Boisse et al., 1998)) may 
render the gaseous medium in damped Lyman-a systems 
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only more extreme versions of those in local dwarf systems 
like the LMC and SMC ^ But it has also been argued that 
the apparent recognizability of patterns in the absorption 
spectra has been over-interpreted (Izotov et al., 2001), and 
that the similarity of intermediate and low ion kinematics 
(Al III and C II or Fe II; C IV and C II behave differently) 
could mean that substantial ionization corrections to the 
metallicity are needed (however, see Vladilo et al. (2001) 
for a contrary opinion). 

At the present time there seem to be several lines of 
evidence suggesting that, unlike the local ISM (where low- 
altitude neutral gas is perhaps 2/3 cool and 1/3 warm 
and the overall ratio including high-altitude material is 
1/2 and 1/2) the gas in damped Lyman-a systems is 
more predominantly warm. Examples include the high H 
I spin temperatures inferred from comparison of A21cm 
and Lyman-a absorption (Wolfe & Davis, 1979; Carilli 
et al., 1996a, c; Chengalur & Kanekar, 2000; Kanekar & 
Chengalur, 2001) (but see Lane et al. (2000)), the similar- 
ity of low and intermediate ion kinematics mentioned in 
the preceding paragraph, and the very low column den- 
sities of II2. Along lines of sight with reddening F,b-v> 
0.05 mag (N(II) > 3 x 10^*^ cm~^) in local Copernicus 
spectra, it is always the case that a few percent or more of 
the neutral hydrogen is molecular. By contrast, the molec- 
ular hydrogen fraction in damped Lyman-a systems is 2-4 
orders of magnitude lower (i.e.lO"'* — 10~®), which can be 
interpreted as meaning that the gas temperature must be 
above 3000 K (Petitjean et al., 2000). 

Here we consider the formation of molecular hydro- 
gen in a gas which is in two-phase thermal equilibrium 
at low metallicity. Perhaps because II2 has been observ- 
able so rarely in the local ISM there is not a big literature 
on this subject, but the extant 1970's-era Copernicus ob- 
servations are well-explained in this way (Liszt & Lucas, 
2000) using modern shielding factors for radiative dissoci- 
ation (Lee et al., 1996). The only surprise (if there indeed 
is one) is the low densities that are required to start II2 
formation locally, and the fact that even a "standard" H 
I cloud (Spitzer, 1978) should have a molecular fraction 
of 10-30% deep inside. Local diffuse clouds also have sur- 
prisingly high abundances of complex polyatomic species 
which follow immediately upon the presence of II2 (Liszt 
& Lucas, 1996; Lucas & Liszt, 1996), a phenomenon which 
is not well understood but which can be used to account 
for the observed abundancess of simpler species such as 
CO (Liszt & Lucas, 2000). 

In the course of this work, we consider the ionization 
and fine-structure excitation of carbon, and the inferences 
which may be drawn from observations of carbon in the 
damped Lyman-a systems. While the importance of car- 
bon to the physical state of the gas and the conditions for 
forming II2 cannot be stressed too highly, the extended 
discussion here in fact arose because of something of a 

^ Boisse et al. (1998) show that there there is a bias against 
finding high metallicity damped Lyman-a systems at high N(H 
I), because the background sources are extinguished. 



coincidence. Searching the literature, it quickly became 
evident that there is a high degree of overlap between the 
two relatively scant datasets for carbon (chiefiy, C I, C II, 
and C II*) and H2. 

The plan of the current discussion is as follows. In Sect. 
2 we lay out the basics of a calculation of two-phase equi- 
librium essentially following Wolfire et al. (1995). We ex- 
amine the sensitivity of two-phase equilibrium to varia- 
tions in abundance, depletion, incident radiation and the 
like, in order to extract from the calculations those as- 
pects which can be related to existing absorption line data 
on damped Lyman-a systems. In Sect. 3 we compare the 
results of these calculations with observed hydrogen and 
carbon column densities: we show that the observations 
are consistent with an origin in largely cool gas for the 
damped Lyman-a systems at lower z (z < 2.3) and in 
warm gas at higher redshift. In Sect. 4 we describe calcu- 
lations of the II2 abundance in both warm and cool diffuse 
gas, under conditions of varying metallicty, etc., employ- 
ing (slow) gas-phase processes to explore the minimum 
expected amounts of molecular gas, and the more usual 
grain-catalysis (as in Liszt & Lucas (2000)) for cooler re- 
gions of higher molecular fraction. We also compute the 
variation of the molecular hydrogen fraction in small gas 
clots of constant density. From this, it follows that the low 
metallicities of damped Lyman-a systems are by them- 
selves sufficient to cause decreases in the molecular frac- 
tion by many orders of magnitude, even if cool neutral 
clouds are present. 

2. Two-phase equilibrium of the ISM 

2.1. Particulars and parameter sensitivities 

The ISM is so complicated that any model of it is bound 
to be heavily idealized and stylized, especially the assump- 
tion of a strict equilibrium. But some aspects, especially 
the presence of phases - discrete regimes of density and 
ionizat, seem robust to variations in the underlying pa- 
rameters. The existence of multiphase gas seems to be 
quite general, highly conserved across space and time. 

The basis of the present work is a calculation of two- 
phase equilibrium following work on the local ISM by 
Wolfire et al. (1995). Notable constraints on the model lo- 
cally are the thermal pressure range P/k = 10^ — lO'* cm~^ 
K observed (via C I and C I*) in neutral gas locally by 
Jenkins et al. (1983), and the electron density in warm 
gas at the Solar radius Ug ~ 0.02 cm~^ (Taylor & Cordes, 

1993) . As noted by Wolfire et al. (1995), strict two-phase 
equilibrium is an idealized goal toward which the ISM may 
tend, but true equilibrium is hard to attain, especially in 
warm gas where the cooling timescales are long. 

Basically, this "standard" model of the local diffuse 
ISM is driven by the heating due (in increasing order of 
importance in warm gas) to cosmic-rays, soft X-rays, and 
- the dominant mechanism in both warm and cool gas - 
the photoelectric effect on small grains (Bakes & Tielens, 

1994) . The charge state of the small grains, their heating 
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Fig. 1. Ionization and thermal equilibrium calculations in atomic gas. The pressure P/k is shown as a function of the density of 
H-nuclei n(H), for four sets of parametric variations, a) Upper left; the metallicity, (dust/gas, C/H, O/H etc.) varies in steps of 
2 from 4 times to 1/32 times its reference value; b) Upper right, carbon and oxygen are removed (depleted) from the gas phase 
in steps of 2; c) Lower left, the "interstellar" (ambient) radiation field (ISRF) varies from 4 times to 1/32 times the reference 
value; d) Lower right, the impinging flux of soft X-rays is scaled. 



rates and the recombination of ions on them (which dom- 
inates over gas-phase processes at moderate densities in 
diffuse gas) are all very sensitive to the electron density, 
so local thermal and ionization equilbria are tightly cou- 
pled. Given that carbon is the main source of cooling in 
cooler gas (O I to a lesser extent) and an important source 
of free electrons, it follows that the equilibrium conditions 
and the condition of carbon in the gas are very tightly 
coupled. 

We recently used this code (Liszt, 2001) to discuss the 
behaviour of the H I spin temperature in local warm gas 
(the H I there is not generally thermalized by collisions in 
multiphase equilbrium) and to calculate the abundances 
of H2 and CO in denser diffuse gas (Liszt & Lucas, 2000). 
The latter reference describes the model in more detail 
than will be given here, especially the various collisional 
processes which drive fine-structure cooling in the gas 
(they are, of course also discussed in the original refer- 
ence by (Wolfire et al., 1995)). 

Our basic calculation differs from that of Wolfire et al. 
(1995) in two minor ways which are subsumed by the ex- 



tent of our parameter variations. The earlier authors took 
the locally-determined soft X-ray flux from Garmire et al. 
(1992) and decomposed it into three components, each 
of which they represented by a (physically-motivated) 
plasma emissivity which was then separately attenutated 
(or not) by an assumed column of warm gas; the three 
contributions were then summed to provide the assumed 
incident X-ray spectrum at a typical location. Their final 
spectrum is characterized by a quantity they called Nw, 
the overall column of warm gas attenuating part of the 
incident spectrum ^ whose standard value they took to be 
Nw = 10^^ H atoms cm . We used a simpler representa- 
tion of the soft X-ray flux, whereby the observed spectrum 
is attenuated by a column of neutral gas Nw directly. Our 
standard value of Nw = 2 x 10^® H atoms cm~^ produces 
a pressure-density curve which differs little from that of 
Wolfire et al. (1995). 



^ The soft X-ray flux is so fragile that any cool gas interposed 
would absorb all of it, given that cool gas does not appear in 
small increments 
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Fig. 2. As in Fig. 1, but the quantities shown in each panel are the energy loss rate nA due to radiation in the C 11* A158jU 
cooling transition (lower, full curves) F, the total energy input to the gas (upper, shaded curves), both in units of erg s~^ H~^. 
In cool gas (to the right), C 11 is by far the dominant coolant; in warm gas, most of the cooling occurs in the Ly-a line. 



We also used a somewhat different set of reference 
(Solar) atomic abundances whereby log(C/H) = -3.44, 
log(0/C) = 0.32, namely those which accompanied the 
distribution of the soft X-ray absorption cross- sections 
from (Balucinska-Church & McCammon, 1992), which 
both we and Wolfire et al. (1995) used (in their most re- 
cently updated version). In both calculations the reference 
model has (following Wolfire et al. (1995)) no depletion of 
the gas phase oxygen and carbon onto grains, which might 
be acceptable for warm gas (where carbon does not bear 
the brunt of the cooling or contribute many of the elec- 
trons, see below) but is unlikely to be reasonable in local 
cool gas; if there are grains, they have to be made of some- 
thing! The gas phase depletion is considered an adjustable 
parameter here, independent of the metallicity. 

The net effect of these two differences is slight. With 
the exception of one system, the observations show that 
carbon is present in the gas in the same proportion as 
other species which are typically undepleted locally, i.e. 
there is low metallicity, but little global carbon depletion. 

In Fig. 1 we show the basics of two-phase equilibrium 
and its sensitivity to some assumed parameter variations. 
In these curves of pressure vs. density (one assumes a 



density, calculates the ionization and thermal equilibrium, 
sums the particle densities and multiplies by the derived 
temperature), the regions where dP/dn < are held to be 
unstable and unlikely to occur in nature. Those pressures 
for which there are two densities n' having dP/dn|„/ > 
are those for which "two-phase" equilibrium is possible. 
The curve labelled "1" in all panels is the same (reference) 
model; perhaps most clearly at upper right it is apparent 
that the reference model provides for two-phase equilib- 
rium over precisely the pressure range which is observed 
in local gas (Jenkins et al., 1983). Typically (but see be- 
low and the lower-left panel), only warm (7000 K - 10000 
K) neutral gas appears if the density and pressure are be- 
low the two-phase regime and only cold (below 1000 K) 
neutral gas appears if the pressure and density are larger. 
If conditions are such that the warm gas may persist up 
to higher density and pressure, or when the cold medium 
cannot exist except at higher density and pressure, the 
medium is more likely to be warm. 

At the top left, the metallicity is varied by factors of 
4 to 1/32 relative to the reference value. Many properties 
of the model actually vary in concert including; the ele- 
mental and gas-phase C and O abundances (there is no 
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depletion); the columns of all other species heavier than 
Helium providing soft X-ray absorption; and the number 
of grains both large and small (the latter scales the heating 
rate of photoelectric ejection). Vladilo (1998) shows that 
the dust/metal ratio D/Z « 0.6 at higher redshift rela- 
tive to conditions in the Milky Way (MW) so that D/G = 
D/Z X Z/G changes mostly because of metallicity (Z/G), 
not D/Z. In fact, such an overall scaling has a surprisingly 
modest effect compared to some other changes discussed 
next. The warm phase may persist up to slightly higher 
pressures and densities in gas of lower overall metallicity 
but the existence of two-phase equilibrium is not threat- 
ened as long as the balance between heating (small grains) 
and cooling (carbon) is maintained in denser gas (as it is 
in this case). 

At upper right, the gas-phase abundances of C and O 
are decreased while all other parameters are held fixed and 
this actually has a much more profound effect than scal- 
ing the grain (heating) and gas (cooling) abundances in 
tandem. Carbon provides some heating of the gas (via its 
ionization, see Spitzer & Scott (1969)) but is more impor- 
tant as a source of electrons, infiuencing the grain heating 
and recombination rates, and as a coolant. Clearly, when 



the ambient fiuxes are maintained while the main coolants 
(carbon and oxygen) and source of electrons are depleted, 
it is substantially more difficult to provide enough pres- 
sure to maintain gas in the cool phase. This hearkens back 
to our discussion, slightly earlier, of the lack of inclusion 
of depletion in the reference model for MW conditions. In 
order to maintain two-phase stability over the reference 
range of pressures in the presence of gas-phase depletion, 
the incident optical/uv ionizing radiation field must be de- 
creased somewhat. In fact this might not be unreasonable 
locally for cool gas, which does not occur in arbitrarily 
small increments and therefore exists under conditions of 
non-trivial extinction (especially at uv wavelengths). At 
lower metallicity, this sort of compensation is less obvious 
and (consistent with some aspects of the observations) the 
pressure may in fact be somewhat higher in cool gas. 

At lower left, we scale the ambient optical/uv radiation 
field which ionizes carbon (directly) and provides heating 
and electrons via the photoelectric effect on small grains. 
Regions whose ambient radiation field (ISRF) falls even a 
factor four below the average will be predominantly cool. 
This is an inversion of the normal order of things in diffuse 
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4. Fine-stmcture excitation calculations for varying 
z'x2.73 K. Top: the ratio of first-excited to ground- 



Fig 

Tcmb 

state populations in neutral carbon, CI. Middle, the ratio of 
excited to ground state populations in once-ionized carbon, C 
11. Bottom, as in the middle panel but for the case where car- 
bon is depleted by a factor 4 from the gas-phase, relative to 
the standard model. 



gas and darker clouds in the nearby ISM occur at higher, 
not lower, thermal pressure. 

Last, at lower right, we scale the soft X-ray flux. Cool 
regions of higher density will probably never see the full 
extent of the ambient soft x-ray flux because even a small 
region of appreciable density accumulates a high degree 
of X-ray absorption. The heating due to soft X-rays tends 
to increase on a per-event basis as the X-ray spectrum 
hardens somewhat after absorption, and tends to decrease 



— as the electrons cause more secondary ionization and less 

- direct heating when the ionization fraction is smaller (in 

- neutral gas). 

2.2. Thermodynamics 



Fig. 2 shows, for the same set of parameter variations, the 
total gas heating rate F and the loss of energy (nA) due 
to cooling in the A158// C II line (discussed below at some 
length), in units of ergs s~^ H"^. For the reference model, 
some 90% of the heating is due to the photoelectron ef- 
fect operating on small grains, see Fig. 3 of Wolfire et al. 
(1995). In warm gas, the burden of cooling is taken up 
by excitation of the Ly-a lines under the highly-idealized 
assumption that all photons escape; in fact, these pho- 
tons may travel substantial distances but by and large do 
not actually escape a two-phase galactic layer before be- 
ing absorbed by dust, at Solar metallicity: the situation 
may actually be very different in three-phase models, see 
Neufeld (1991) and Liszt (2001). In cool gas, the cooling 
is due almost entirely to carbon and oxygen fine-structure 
excitation, in the reference model. The curves of nA in 
Fig. 2 directly show the brightnesses (normalized; per H) 
of the C II* A158// and O I* cooling transitions. 

There is really a quite profound change in the lo- 
cal thermodynamics when the metallicity or ionizing fiux 
varies (Fig. 2, upper and lower left). In particular, while 
the energy input per H into cool and warm gas is nearly 
the same for the "Solar" metallicity and the standard 
ISRF or higher, it is much smaller in cooler gas when 
the metallicity is low. In the Milky Way, the carbon cool- 
ing rate can be used to infer the heating in both cool and 
warm gas. In systems of low metallicity such is not the 
case. 

In general, the strength of the C II A158// line varies 
with metallicity at all densities; with carbon depletion at 
most but especially at all lower densities; with the strength 
of the ISRF in cool gas; and only slightly at low density 
with variations in the soft X-ray fiux. From the results at 
upper right we see that C II may not be the dominant 
source of cooling in cooler gas, at quite high densities, 
if the depletion (not the metallicity) is extreme; oxygen 
becomes the preferred coolant due to the higher energy 
separation in its ground state fine-structure levels. There 
is no reason to believe that such conditions occur widely in 
diffuse gas, either locally or in damped Lyman-a systems, 
although one source in Table 2 (PKS0250) seems quite 
deficient in carbon, given its quoted metallicity. 



2.3. Carbon ionization 

Fig. 3 shows the ratio of the two lowest ionization states 
of carbon n(C II)/n(C I). It differs by a factor of 40 - 
100 between warm and cool neutral gas for the standard 
model, and is relatively insensitive to parameter changes, 
redshift, etc. This makes it a sensitive indicator of the 
thermodynamic state of the g discussed in Sect. 3, 
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but also causes confusion when mixtures of the two phases 
are observed along the same line of sight (see Sect 2.6). 
Some aspects of the behaviour shown in Fig. 3 are counter- 
intuitive; for instance, a stronger ISRF maintains a high 
C II/C I ratio into denser gas, but the C II/C I ratio is 
actually smaller in very tenuous gas when the radiation 
field is higher. This arises in part because the equilibrium 
pressure increases with the strength of the radiation field 
and partly due to the level of ionization in the gas (which 
is not solely determined by carbon). 

If the C II lines are saturated or N(C II) is other- 
wise unknown, but N(C I) and N(II) are available, the 
thermodynamic state of the gas may also be inferred by 
computing N(C) /N(H) ^ N(C II) /N(H) N(C II) /N(C I) 
X N(C I)/N(II). Presumably, if the gas is actually cool, 
the inferred carbon abundance will be much too large (in- 
consistent with the metallicity) if the high n(C II)/n(C I) 
ratios typical of warm gas are assumed. 

2.4. CI fine-structure excitation 

The triplet fine-structure of the C I ground state is also 
accessible to spectroscopy. It is actually rather insensitive 
to the parameter variations noted above (with the slight 
exception that C I*/C I increases modestly with increasing 
X-ray fiux at low density) but responds very strongly to 
a change in the cosmic background level because the first 
excited state is only 23.2 K above ground. Fig. 4 (top) 
shows the ratio of populations in the lowest two levels as a 
function of Tcmb- Clearly, this ratio is a sensitive indicator 
of the ionization state at low redshift (it differs by a factor 
of about thirty between warm and cool gas) but is rather 
insensitive for redshifts above 1.5 or so. The C I*/C I 
ratio measures Tcmb at high redshift in the context of 
our models (see Roth & Meyer (1992)). Although perhaps 
within observational errors, it follows from the results in 
Fig. 4 that quoted column densities of neutral carbon for 
high-z systems should include a correction factor for the 
excited state population. 

2.5. Cooling rates and C II fine-structure excitation 

The rate of energy loss due to a radiative transition be- 
tween two levels k and j having statistical weights gk and 
gj , level populations n^ and nj (units of cm~^), energy 
separation E^j and spontaneous emission rate A^j (s""'^), 
immersed in an isotropic radiation field characterized by 
Tcmb is (Goldreich & Kwan, 1974) 

riAkj = ylkj/3kj-E'kj(«k - ("jgk/gj - nk)Fkj (Tcmb)) 

where Fkj (x) = l/(exp(Ekj/kx)-l). /3kj is the probabil- 
ity of photon escape, typically varying as /3kj = (l-exp(- 
urkj))/(urkj) where u is a geometrical factor of order unity 
and Tkj is the optical depth of the transition. The escape 
probabability is taken as unity in these two-phase mod- 
els, which is acceptable because the various fine-structure 



transitions do not achieve high optical depth in typical 
diffuse clouds. In calculation of the local statistical equi- 
librium, these expresssions allow for a proper evaluation 
of the fine-structure cooling rates for varying Tcmb • 

The 133.5 and 133.6 nm C II and C II* absorption 
lines have received more attention because of the obvious 
possibility that C II is the dominant ionization stage. But 
discussions of the excitation of the C II* fine-structure 
levels often fail to recognize the role of carbon in deter- 
mining the conditions under which it is observed; in many 
cases (see Fig. 2) C II will be the dominant coolant when 
it is the dominant ionization stage. This means that the 
brightness of the C II A158// line and the C II*/C II ra- 
tio are determined by thermodynamic equilibrium - the 
energy input to the gas which occurs at a given density 
and temperature. So, if the energy input into the gas does 
not change, neither will the brightness of the A158// line 
or the strength of the 133.5 nm C II* absorption line. In 
Fig. 2 at upper right, for a density of n(II) = 20 cm~^, 
it requires a factor 8 drop in the the amount of carbon in 
the gas to weaken the A158// line by about a factor 2, all 
other things being equal. 

Fig. 5 shows the parameter sensitivities of the C II*/C 
II ratio. Weak variation in the n(C II*)/n(C II) ratio in 
cool gas is shown in many ways in the panels of Fig. 5. For 
the reference model, the ratio varies by only 50% between 
n(H) = 2 cm-3 (Tk= 1000 K) and n(H) 40 cm'^ (Tk= 75 
K). The variation with carbon depletion at upper right il- 
lustrates the effects of thermodynamic considerations. The 
ratio n(C II*)/n(C II) increases as carbon is removed from 
the gas or n(C II) decreases, because the brightness of the 
A158// line (and n(C II*)) must stay constant in order to 
carry away a given amount of energy which is being in- 
put to the gas (see Fig. 2). The behaviour at lower left, 
where the C II*/C II ratio is a scaled representation of 
the cooling rate (also see Fig. 2) shows how the brightness 
of the A158// line, or the amount of excited-state C II*, 
must change proportionally as the energy input to the gas 
changes. 

In warm gas where its excitation is weak and C II is 
not the dominant coolant, the population of the C II* 
level is susceptible to change with Tcmb - This is shown in 
the lower two panels of Fig. 4. The ability of the C II* 
population to discriminate between phases may be lost at 
higher redshift. 

In this context it is important to note that the AE/k 
= 91 K energy separation of the C II fine-structure levels 
(to a lesser extent the 232 K separation in O I), is ab- 
solutely crucial to the similarity of two-phase equilibrium 
conditions out to moderate redshift, say z < 10. If the 
gas locally were cooled by softer photons of energy AE 
such that we could observe out to some redshift where 
AE/kTcmb^ 4, the distant gas would have to have a dif- 
ferent coolant. It is apparent from the behaviour of the C 
I*/C I ratio in Fig 4 that C I*, for example, could not be 
an important gas coolant at z > 3. The same is true of 
the lower lines of CO (where AE/k = 5.5 K, 11.0 K, 16.5 
K, etc.). 
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Table 1. Hydrogen and carbon at lower redsliift 



Source 


3C286 


0454+039 


1756-1-237 


1331-1-170 


0013-004 


1232-1-0815 


PHL957 


z 


0.69 


0.86 


1.67 


1.78 


1.97 


2.24 


2.31 


[Zn/H] 


-1.22 


-1.1 


(-0.8)" 


-1.27 


-0.80 


-1.20 


-1.38 


N(H 1) 


1.8E21 


4.2E20 


2.0E20 


1.5E21 


5.0E20 


7.9E20 


2.8E21 


N(H2) 










6.9E19 


6.3E16 


<5.0E15 


N(C 1) 


>4E13 


4.4E13 


1.4E13 


1.7E13 


6.3E13 


4.0E13 




N(C 11) 








4.5E16 


2.7E16 




7.8E15 


N(C 11*) 

N(C ll)n,'' 






7.9E13 


1.4E14 


1.9E14 




7.6E13 


3.8E16 


1.2E16 


(1.1E16) 


2.9E16 


2.8E16 


1.8E16 


4.1E16 


N(C 11)/N(C 1)^= 


(<952) 


(293) 


(803) 


2647 


446 


445 


>1734 


N(C 11*)/N(C 11) 






(0.0072) 


0.0031 


0.0070 




0.0097 



" assumed for illustrative purposes; other (entries) follow 

N(C ll)m is the amount of carbon if [C/Zn] is Solar 
^ values in parentheses use N(C ll)m 
References: 



3C286: Boisse et al. (1998); Roth & Meyer (1992) 
0454-1-039: Boisse et al. (1998); Steidel et al. (1995) 
1756-1-237: Roth & Bauer (1999) 
1331-1-170: Chaffee et al. (1988); Ge et al. (1997) 
0013-004: Ge et al. (2001) 

1232-1-0815: Black et al. (1987); Srianand et al. (2000) quote N(H2) = 1.5E17 cm"^ 
PHL957: Black et al. (1987) 
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Fig. 6. Molecular fractions 2n(H2)/n(H) in free space resulting 
solely from slow gas-phase H2 -formation mechanisms, for the 
two-phase models and paramter variations shown at left in Fig. 
1. See Sect. 4 of the text. 

2.6. Observing admixtures of the two phases 

When a property like the N(C II*)/N(C II) ratio - call it 
R - takes on values Rw and Rc in warm and cool gas, 
respectively, and a transparent superposition of phases is 
observed to have a global value Robs , Robs is related to 
the proportion of gas in the two phases as 



Rc 



Rw 



Rw — Rc Robs 
where Fc is defined as the fraction of the gas in the 
cool phase. 

In stable two-phase gas (Fig. 3 at upper left) Rw = 
20000 and (roughly) Rc = 700/n(II) in gas of Solar metal- 
licity (Rc varies rather less at low metallicity). In this 
case, measuring a global average of, say, Robs = 1000 im- 
plies only that there is some fraction Fc< 1 of gas having 
Rc< 1000 along the line of sight. 

This ambiguity cannot generally be resolved with re- 
course to comparison of profiles, which are largely undif- 
ferentiated in the C II line, having been broadened due to 
some combination of; underlying kinematics in the host 
system; doppler-broadening of individual gas complexes; 
and, saturation. In the Milky Way, C II lines are typically 
fit with much larger b- values (Doppler broadening param- 



eters) even when a single cloud or gas complex dominates 
the line of sight to a nearby star (c/. Morton (1975)). 
Because of this one cannot, for instance, try to compare 
the C II* with only that fraction of the C II which overlaps 
it in velocity. 



3. Warm or cool gas in damped Lyman-a 
systems? - observations of C I and C II 

In Tables 1 and 2 we have gathered information on the 
neutral and ionized carbon lines available in the litera- 
ture. For each source we show the observed column den- 
sities, and the metallicity derived from zinc, which is be- 
lieved to be (at most) only lightly depleted (Pettini et al., 
1994; Vladilo et al., 2000). From the metallicity we derive 
a quantity N(C II)m which is meant to represent the max- 
imum amount of carbon available in the gas phase, i.e. it 
is the column density of carbon assuming that carbon and 
zinc have the same metallicity and depletion. Comparison 
with the measured C II column densities shows that car- 
bon is not lacking, with the notable exception of PIIL957. 
N(C II)m is used as a surrogate where measured values of 
N(C II) are unavailable at low redshift in Table 1. 

The ratio N(C II) /N(C I) is much larger for the sources 
at high z in Table 3, above 13,000-30,000, consistent only 
with an origin in overwhelmingly warm gas: the fraction of 
cool gas must be of order a few percent at most. At lower 
redshift, two of the measured ratios, and that inferred for 
the source 1756-1-237, are below 1000, consistent with cool 
neutral gas in the stable region of the two-phase equilib- 
rium, at thermal pressures P/k = 3,000 - 4000 cm~^ K, 
and densities n(H) = 10-70 cm'^. The data for PHL957 
are also consistent with stable cool gas at n(II) ^ 10 cm~^, 
p/k 2500 cm~^ K, while 1331-1-170 is in the marginally 
unstable region at n(II) = 10 cm~^, p/k = 2800 cm~^ 
K. For the systems at lower redshift, there must be a very 
large proportion of cool gas, though we cannot say exactly 
how much. 

The N(C II*)/N(C II) ratios are noticeably smaller at 
high redshift in Table 2, and would presumably be smaller 
still by a factor 2 or more, were it not for the cosmic 
background radiation, which has a noticeable effect (only) 
above z = 2.5 in warm gas (see Fig. 4). The N(C II*)/N(C 
II) ratios observed at lower redshift, in presumably cool 
gas, are 50-100% higher than in the models. This could 
be explained by an enhanced ISRF (Fig. 5, bottom left 
panel) or by some depletion of carbon for PIIL957, where 
the metallicity of carbon is (uniquely) very low compared 
to that of Zn, and N(C II*)/N(C II) is rather large. 

To summarize, comparison of the tabulated observa- 
tional results with the calculated properties of carbon in 
two-phase media shows quite unambiguously that the ob- 
served gas has a substantial contribution from cool gas at 
lower redshift but is very largely warm at z > 2.8. Whether 
this shift from cool to warm gas is systematic or coinci- 
dental remains to be explored. 
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Table 2. Hydrogen and Carbon at highest redshift 



Source 


1337+113 


0528-250 


0347-382 


0000-262 


z 


2.80 


2.81 


3.03 


3.39 


[Zn/H] 


-1.00 


-0.91 


-1.23 


-2.07 


N(H I) 


8.0E20 


2.2E21 


2.52E20 


2.6E21 


N(H2) 


<5.0E16 


6.0E16 


8.2E14 


1.1E14 


N(C I) 


<1.6E13 


<5.9E12 


<4.0E11 




N(C II) 


2.0E17 


1.7E17 


5.1E15 




N(C II*) 




3.6E14 


1.9E13 




N(C II)n, 


2.8E17 


9.6E16 


3.1E15 


7.9E15 


C II/C I 


>12658 


>28862 


>12700 




C II*/C II 




0.00211 


0.00389 




References: 



1337-1-113: Lanzetta et al. (1989) 

0528-250: Ge et al. (1997); Srianand & Petitjean (1998) 
0347-382: Levshakov et al. (2002) 

0000-262: Prochaska & Wolfe (1999); Levshakov et al. (2000) 



4. The formation of molecular hydrogen 

Levshakov et al. (2002) provide a summary of the molec- 
ular fractions observed toward 14 sources (their Table 2) 
including most of those discussed here. The lowest molecu- 
lar fraction is seen toward 0000-262 at z = 3.39, the high- 
est toward 0013-004 at z = 1.97 (see Tables 1-2 here). 
The latter is the only direction for which the molecular 
fraction exceeds 1/3000, and most are below 10~®. In this 
Section we address the existence of damped Lyman-a sys- 
tems with large hydrogen column densities, substantial 
amounts of cool gas, but very small fractions of molecular 
hydrogen. The models of H2-formation presented here are 
very similar to those of Liszt & Lucas (2000), who consid- 
ered the formation of CO in local diffuse gas, given certain 
other conditions like the presence of HCO"*" , but there are 
a few differences which we now remark. 

4.1. H2-formaf;on in the gas phase 

The present calculations include X-ray heating and ioniza- 
tion, because they are integral to the question of two-phase 
equilibrium. X-rays have little effect on the large neutral 
gas columns which harbor appreciable molecular column 
densities nearby in the Milky Way but are also included 
now because we are interested in understanding the mini- 
mum molecular fractions which can be expected, and these 
are presumably set by the slow gas-phase processes (not 
involving grains) which formed the first stars, and which 
occur in low-density regions of (unshielded) free space. 

Gas phase H2-formation occurs via the exothermic re- 
action pathways H -|- e~ — H~ -|- hi>, H~ -|- H — H2 -|- 
e- and H+ H ^ H+ /li^, H+ H ^ H2 H+ . Most 
of the basic reactions are listed by Puy et al. (1993) (all 
can straightforwardly be located in the UMIST reaction 
database) but their discussion of early-universe conditions 
must be augmented by a value for the photodissociation 
of H~ , by inclusion of the endothermic (22000 K) destruc- 
tion of H2 by H"'", by cosmic-ray ionization (here taken as 



C = 2 x~^^ s~^ per H and 1.08 ( per H2), and by the 
photodissocation rate of H2 in free space (which we take 
from Lee et al. (1996)). Current values of the reaction 
rates for all important processes are given in the UMIST 
reaction database, whose values we employed unless oth- 
erwise noted. Species to be followed during modelling of 
the chemistry include H I, He"*", H"*" and e~ - all given by 
the calculations of two-phase equilibrium -as well as H" , 
H2 and H+ . 

In order to formulate a treatment of the variation of 
the photodissociation of H" with extinction, we integrated 
the cross-sections of Wishart (1979) over the local ISRF, 
finding an unshielded rate of 1.5 x 10~^s~^, just over half 
of which arises at wavelengths beyond 800 nm; the photo- 
dissociation is, therefore, not strongly attenuated under 
the conditions discussed here and we elected to ignore ex- 
tinction in this regard. The rate quoted in the UMIST 
database is 2.4 x 10~^ s""'^, which we used in our calcula- 
tions. 

Fig. 6 shows the free-space abundance of II2 arising 
solely from gas-phase processes in the two-phase models 
under conditions of varying metallicity and ISRF (the two 
left-hand panels of Fig. 1). Unlike grain formation scenar- 
ios, which take advantage of high N(II) to boost the molec- 
ular fraction at high density, free-space gas-phase forma- 
tion of II2 does not seem to distinguish between warm 
and cool conditions. The molecular fractions calculated in 
Fig. 6 correspond well with the smallest values in the lo- 
cal ISM, or, for that matter, in damped Lyman-a systems. 
Unfortunately, this minimum is not diagnostic of the host 
gas conditions. 

4.2. The H2 abundance in cool clouds 

II2 formation in cool neutral gas clouds is illustrated in 
Fig 7, which indeed shows why even damped Lyman-a 
systems with appreciable cool gas still may lack molecu- 
lar hydrogen. To create this diagram, we considered (fol- 
lowing Liszt & Lucas (2000)) a spherical clot of gas of 
constant density, immersed in isotropic radiation fields 
(X-ray, cosmic-ray, optical/uv, etc.). This was compu- 
tationally divided up into 128 radial zones, in each of 
which we derived the temperature/ionization structure 
and II2 abundance. The latter requires iteration, because 
the maintenance of II2 is a sharply non-linear process de- 
pendent on the II2 column and extinction between any 
point and free-space (Lee et al., 1996). We adopted a fairly 
straightforward relaxation method which converged with 
gratifying rapidity. 

Calculation of the abundance of molecular hydrogen 
is typically made feasible by employing a set of shielding 
factors which account in an average way for the many very 
complicated effects of line-overlap and radiation transport 
in the dissociation process. We used the shielding factors 
of (Lee et al., 1996) which were calculated for local gas. 
The justification for this is that the dominant effect requir- 
ing consideration here is the order of magnitude change in 
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Fig. 7. Molecular fraction as a function of fractional distance 
into cool clouds of density n(H) = 32 cm~^ and various central 
column density N(H) as indicated. The uppermost curve in 
each panel is for local conditions. The labelled curves are for 
varying metallicity and a dust/metal ratio D/Z = 0.6 (60% of 
the local value). The bottom curve represents formation solely 
in the gas phase, in the limit of zero metallicity parameter; see 
Sect. 5 of the text. 



the number of grains at very low metallicity, not the factor 
of two differences in individual grain properties between 
local grains and those seen, for instance in the Magellanic 
clouds. The parametrization of Pei (1992) for the Milky 
Way, LMC and SMC shows that, for a given amount of 
B-band extinction, the grain distribution provides suc- 
cessively somewhat more extinction at (say) Ly-a as the 
metallicity declines; the inference is that graphite grains 



disappear and silicates do not. But this effect is dominated 
by the overall diminution of the extinction with lowered 
metallicity. 

Figure 7 shows the radial variation of the fraction of 
H-nuclei in molecular form over spherical gas clots of con- 
stant density n(H) = 32 cm~^ for different column densi- 
ties N(H) through the center of the clot. The mean line of 
sight averaged over the circular face of such a body inter- 
sects it at an impact parameter of 1/3 of the radius (at a 
value 0.33 along the horizontal axis in Fig. 7), where the 
column density is 3/4 of that through the center. 

In each panel of the Figure there are 8 vertically- 
separated curves. At top, shaded, is the result which would 
apply in the Milky Way, where we have taken the dust/gas 
ratio as observed locally (the reference model of the two- 
phase calculations) and depleted carbon and oxygen in 
the gas phase by a factor of 2.4. The bottom curve, also 
shaded, is the result when the metallicity goes to zero and 
only the gas-phase formation of H2 is included; a modest 
amount of self-shielding occurs and the molecular fraction 
is slightly higher than in free space (Fig. 6). The interme- 
diate curves assume a dust /gas ratio 0.6 relative to the ref- 
erence model for the Milky Way following Vladilo (1998) 
(a small effect at higher N(H) but of real importance to 
the thinnest model) and no depletion of carbon and oxy- 
gen. These curves are labelled by their varying metallicity 
as in the previous diagrams. 

The cloud with N(H) = 4 x 10^° cm~^ would be a com- 
pact (4 pc), cool (130 K) Spitzer (1978) "standard" cloud 
in the Solar neighborhood. It would also be very substan- 
tially molecular if found in the Solar vicinity, but a factor 
4 decline in the metallicity suffices to reduce the molecular 
fraction by some four orders of magnitude. Even the model 
having a four times higher column density (compare with 
the entries in Tables 1-2) cannot sustain an appreciable 
molecular fraction when the metallicity is reduced by a 
factor 10, which is hardly extreme for one of the damped 
Lyman-a systems. 

4.3. The role of geometry 

The role of geometry can also be inferred from Fig. 7. At 
any given metallicity, a cloud with lower N(II) produces 
much less than one-fourth as much II2 as that illustrated 
in the next-lower panel. This is another reason why the 
molecular fraction may vary widely between two lines of 
sight with similar N(II), N(C II)/N(C I), and/or metallic- 
ity (for example). Molecular hydrogen will readily popu- 
late a region when the circumstances are propitious, but 
can easily be prevented from forming by the vagaries of 
local source structure. 



4.4. PKS0528-250 

This source (Table 2) has an overall molecular fraction 
5.5 X 10~® despite the lack of evidence (in carbon) for 
any appreciable amounts of cool gas; earlier we asserted 
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that (very roughly) no more than a few percent of the 
gas could be cool. Carilli et al. (1996b) did not detect 
21cm H I absorption, placing a 2 — a upper limit N(H I) 
< 2 X 10^^ cm~^Tsp. So, the molecule-bearing gas must 
be cool, occupying roughly 1% of the total gas column for 
Tgp ^ 100 K. In the context of our models the gas must 
also be fairly dense, n(H) ^ 100 cm~^, occurring over only 
a very small fraction of the path length (1 kpc or more ) 
occupied by the g whole. 

5. Summary 

The thermodynamics of the ISM near the Sun are strongly 
influenced by the disposition of carbon, in grains and in 
the gas. Spectroscopy of neutral and ionized carbon affords 
the opportunity to probe the processes which are most ba- 
sic to the structure of the gaseous medium, nearby and in 
damped Lyman-a systems. For whatever reason, the very 
distant gas is easily understood in the same terms as that 
seen nearby. Here, we have shown that one seemingly dis- 
parate aspect of high-z systems, their small fractions of 
molecular gas, can also be easily understood. Even when 
cool gas is present, which must be the case for the systems 
discussed here at z < 2.3, abundances of II2 are suppressed 
by many orders of magnitude at lower metallicity as a re- 
sult of the sharply non-linear nature of the processes re- 
quired to maintain substantial columns of II2. No whole- 
sale reorganization of the gaseous medium need be hy- 
pothesized to account for low II2 abundances. Conversely, 
we showed that the slow gas-phase processes which formed 
II2 in the early Universe provide for a minimum molecular 
fraction in the range 10~® — 10~^. 
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